Topological Wigner Crystal of Half-Solitons in a Spinor BEC 
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We consider a one-dimensional gas of half-solitons in a spinor Bose-Einstein condensate. We 
calculate the topological interaction potential between the half-solitons. Using a kinetic equation of 
the Vlasov-Boltzmann type, we model the coupled dynamics of the interacting solitons. We show 
that the dynamics of the system in the gaseous phase is marginally stable and spontaneously evolves 
towards a Wigner crystal. 



Since the early days of quantum mechanics, it is com- 
monly accepted that the Wigner crystal is one of the most 
simple yet dramatic many-body effects. In the seminal 
work published in 1931 [1 , Wigner showed that as a re- 
sult of the competition between the long-ranged potential 
and kinetic energies, electrons spontaneously form a self- 
organized crystal at low densities, in a state that strongly 
differs from the Fermi gas. Experimental observations of 
this effect have been reported in carbon nano-tubes [2]. 
The phenomenon has been discussed in ultracold Fermi 
gases with dipolar interactions [21 H] and also investigated 
in systems with short-range interactions [5] |H] . The con- 
cept of Wigner crystal is used for the description of more 
exotic systems, such as holons (charge solitons) in the 
vicinity of a metal-Mott insulator transition [7] , and the 
ground-state properties of nuclear matter [5]. 

In this Letter, we study the behavior of a one- 
dimensional gas of dark half-solitons (HSs) . HSs are the 
elementary topological excitations of a one-dimensional 
spinor Bose-Einstein condensate with spin-anisotropic in- 
teractions [9]. We derive the interaction potential be- 
tween HSs. We show that in some density range, the 
uncorrelated HS gas becomes unstable and undergoes 
Wigner crystallization, as a result of the competition be- 
tween the fermionic statistics and interactions. Wigner 
crystal is usually associated with long-range interactions, 
such as the Coulomb interaction of electrons. Although 
interactions between HSs are short-ranged, we show that 
their specific form allows the formation of a quasi-ordered 
state exhibiting the quantum signature -density cor- 
relations) of a usual Wigner crystal [10]. We make use 
of a kinetic theory and perform numerical simulations 
to demonstrate the dynamical crystallization of a two- 
species HS gas. We compute the spectrum of the or- 
dered phase and observe that it contains both acoustic 
and optical modes, showing that the optical gap is a func- 
tion of the spin anisotropy. By estimating the density- 
density correlation within the Luttinger liquid theory, we 
demonstrate that the ordered state indeed corresponds to 
a topological version of the Wigner crystal. This ordered 
state could be used to model the low-energy features of 
topological excitations in spinor BECs in more general 
situations. For example, it makes possible the investiga- 
tion of many-body effects in soliton tunneling, which is a 



spectacular yet not obvious quantum phenomenon in the 
nonlinear regime 11-13 . Additionally, there is a growing 
interest in the physics of HSs, mainly due to their spin 
textures, allowing them to behave as magnetic charges 
in the presence of an effective magnetic field P3] • There- 
fore, the study of both the ground-state and dynamical 
properties of an ensemble of HSs can be of a practical 
interest for the future development of magnetricity 15J . 

In what follows, we consider that HSs are traveling 
solutions of the ID spinor Gross-Pitaevskii equation [16] 



m is the boson mass, ai the interaction constant be- 
tween particles having the same spin, and «2 the inter- 
action constant between particles having opposite spins. 
When one soliton is present in a spin component and 
absent in the other, the corresponding object is called 
half soliton. The stability of HSs is ensured in the case 
of spin-anisotropic interactions, when |evi| 3> \a2\ which 
is typically realized in exciton-polariton BEC [TTJ, ITS] . 
Let us consider a HS formed in the er + spin projection, 
for definitencss. In that case, the single-soliton solution 
traveling with speed v is given by 
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where /3 — v/c s , c s ~ \J a\m/c s is the sound velocity and 
£ = h/^/2mainQ denotes the healing length. From me- 
chanical arguments (see the supplemental material [19] 
for details), we can derive the pseudo-potential associ- 
ated with the interaction between the two a+ solitons 
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where represents the effective mass of the solitons. 
In fact, it reduces to an effective potential for the case of 
almost black (v ~ 0) collisions, which reads 
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FIG. 1: (color online) Top: dispersion relation of a HS 
gas exhibiting dynamical instability, obtained for 7+ = 0.1, 
7_ = 0.08 and |Aj = 0.2. The black (purple) full lines cor- 
responds to ui- modes in Eq. pj. The dashed lines 
illustrate to the low-wavelength limit u>± ~ v±k. Bottom: 
phase diagram in the (74., 7_) plane for A = 0.2, depicting 
the gaseous, single-mode and two-mode crystalline phases dis- 
cussed in the text. 



In order to incorporate the interaction with the er_ spin 
projection, we have to take into account the spinor na- 
ture of the ground state of the condensate. In that case, 
the function of two HSs in different components located 
at positions x\ and X2, traveling with opposite speeds v 
and — v can be written using the center-of-mass and rel- 
ative coordinates £ = (x\ + x-i)j1 and r\ = {x\ — X2)/2 
respectively [35]. Integrating over the soliton centroid 
£ to compute the energy E — J £ d£, with £ denoting 
the energy density (see [TH] for details) , we can obtain the 
following potential between almost black 0+ — cr_ solitons 
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where |A| = |aa|/ai is the measure of the spin-anisotropy 
of the interactions. This potential is repulsive for attrac- 
tive inter-spin interactions (012 < 0), which is the case 
considered here. Qualitatively, each HS creates a density 
dip in the other component, and this dip repels the other 
HS because of the negative mass of the latter. 

To describe the dynamics of a gas of HSs, we pos- 
tulate that the phase-space distributions f ± (x,v,t) are 
governed by following kinetic equation of the Vlasov type 

dt 
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where the collision integral is neglected assuming that 
only elastic processes are involved in the system. The 



acceleration term can be simply given by = 
— l/M*dU^/dx, where U^(x) is given in terms of the 
topological potentials Q and ([5| 
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Equation ([6]) generalizes the continuity equation pre- 
sented in Ref. [H] by the introduction of the soliton in- 
teraction potential. We now linearize the system around 
its equilibrium, f ± = Jq + Sf ± to obtain the following 
dispersion relation [19] for the excitations of the soliton 
gas 
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where the integrals I v and Iy have the form 



7±(fc, W ) = ^ ) ^ ^^dv, (9) 
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with X{k) representing the Fourier transform of the po- 
tentials U and V. It is known that the excitations in 
one-dimensional Bose gases below the critical velocity c s 
follow a fermionic statistics, as established in the famous 
Lieb-Liniger theory [25] [26] . Actually, this result is quite 
easy to understand by simply looking at the phase of the 
two-soliton wave function (1) of [19) : exchanging two soli- 
tons located at positions x± and x%, we obtain an overall 
phase shift of 7r in agreement with fermionic statistics. 
Therefore, the equilibrium configuration is that of a one- 
dimensional Fermi gas fo(x, v)^ = Nq /2vpQ(vp — \v\), 
where Vp — nhN^/M* is the one-dimensional Fermi ve- 
locity and Nj^ is the soliton density. In that case, the 
excitation spectrum above the equilibrium configuration 
contains two branches 0J±(k) displaying acoustic behav- 
ior in the long wavelength limit k — > as uj± rs v±k, 
where the velocities are given by cumbersome expres- 
sions of A^ [TS]. In the special case of a symmetric gas 
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with 7 = Ao£ representing the dimensionless soliton con- 
centration, an analogue of the Wigner-Seitz parameter. 
We have here used the relation between the Fermi ve- 
locity of the soliton gas and the sound velocity of the 
condensate, vp/c s — \Z2ttj, which has the meaning of a 
dimensionless interaction parameter. 

Due to the competition between the statistical pres- 
sure and the interactions between the solitons, the dis- 
persion relation Q encodes very interesting features, as 
illustrated in Fig. ([!]). If the soliton gas is dilute enough 
so that the concentration parameter lies below the crit- 
ical value 7* ~ 0.07, the system does not exhibit any 
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ordering, corresponding to a gaseous state. Instead, for 
7 > 7*, the solitons start to perform periodic oscilla- 
tions up to a critical value of the wavevector k* = k*{~f). 
At this point, the frequency softens towards zero and the 
system undergoes crystallization with the lattice constant 
given by d — 2ir/k* . The onset of instability is expected 
at higher k* for higher values of 7, therefore correspond- 
ing to tighter lattices. When the interactions dominate 
and the distance becomes comparable with the charac- 
teristic scale of these interactions, the system forms a 
regular lattice. The ratio of the Fermi velocity to the 
interaction potential ^ provides the following magni- 



tude estimate mf^/Ln 



! I„2 



7 



For the case 



of an asymmetric mixture, a similar analysis allows us 
to conclude about the existence of different phases: a 
gaseous phase (v_ = v+ = 0), defined by the region 
7+ < (16 - |A| 2 - 16%/27r 2 7_)/[167r 2 (V27r 2 7_)], and two 
ordered phases, sustaining single-mode (v- — 0, v + > 0) 
and a two-mode (v-,v+ > 0) oscillations. The latter 
two situations dynamically evolve into a crystalline phase 
sustaining one and two phonon modes, respectively. We 
illustrate the dispersion relation and the onset of insta- 
bility associated to mode softening (Fig. Qa), together 
with the phase diagram (Fig. 0b) , which is zoomed 
in the relevant region of the plane (7+ , 7_ ) where the 
transitions between the different regimes occur. 

In what follows, we analyze the two- mode ordered state 
and show that it indeed corresponds to a crystal exhibit- 
ing the same quantum properties as a Wigner crystal. 
The configuration of minimum potential energy is ob- 
tained for a chain of alternating er + — er_ solitons with 
lattice constant d and quantum oscillations are expected 
to lead only to low-amplitude oscillations around the 
equilibrium distance. For definiteness, we consider the 
half- filling configuration d = 1/Nq. In the harmonic ap- 
proximation, the corresponding Hamiltonian reads 
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where ui (vg) is the deviation of the cr+((7_)-soliton at 
site t from its equilibrium position. The diagonalization 



of (11 1 with the help of Hamilton equations leads to two 
modes, an acoustic mode u)\ and a gapped optical mode 
cj2- In the long- wavelength limit they are given by u>\ &s 
u\k and ui 2 = A 2 + u 2 k 2 , where the velocities u\ and u-i 
are functions of the ratio d/£ — I/7 [19]. Within the first 
neighbor approximation, we obtain 



^ U"(d) V"(d) V"(d) 
u h^d x ——±——, A~2W— — . (12) 



The plasma frequency A is a feature that distinguishes 
the crystal from the Fermi gas, and is a function of the 
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FIG. 2: (color online) Spectrum of a symmetric crystal with 
7+ = 7_ = 0.5, corresponding to the half-filling configuration. 
Acoustic (black) and optic (purple) phonon modes. The inset 
represents the gap frequency A as a function of the parameter 
7- 



spin-anisotropy |A|, A ~ y |A|. In Fig. ([2J), we plot the 
phonon modes wi,2 in the first-neighbor approximation 
for the case of a symmetric crystal. 

In order to investigate the quantum features of the 
crystalline phase, we compute the density correlations. 
Due to the short range character of the interactions, the 
Luttinger liquid theory |27j can accurately describe the 
physical properties of the system. For a model with 



SU(2) spin symmetry such as Hamiltonian ( 11 ) (or, more 
precisely, its second-quantized version) , bosonization pre- 
dicts the decay of the density correlation function as fol- 
lows (251 



K A2 , , A4 
{n x n } ~ + 1+K cos{2k F x) + cos(Ak F x), 

(13) 

where A2 and A4 are some constants. The Luttinger 
parameter [27j [29] K = a/ KfjK v / {K V K V + |A|) can be 
given in terms of the single-spin parameters |19j 



^2^X(ed)[l-cos(2k F £d)], (14) 



with X = U, V. Apart from the x~ 2 dependence, which 
is familiar from the Fermi liquid theory, the properties of 
the system are universally defined in terms of K. There- 
fore, for K > 1/2, the 2kp Friedel-like oscillations domi- 
nate, which is typical of a Luttinger liquid (LL). On the 
contrary, for K < 1/3, 4k F quantum fluctuations dom- 
inate the system, leading to a modulation at the aver- 
age distance d = 1/Nq between the solitons. This corre- 
sponds to a quasi- Wigner crystal (qWC) state, where the 
(quasi) long-range order is due to quantum fluctuations, 
being favored in the low density limit (notice that the 
interaction energy scales as c s /vp ~ 1/7)- In the inter- 
mediary region 1/3 < K < 1/2, there is a mixture (M) 
between the two previous phases. An interesting feature 
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distinguishing qWC stated from the Wigner crystal of 
electrons is the topological nature of the effective poten- 
tials, whose short range is defined by the condensate's 
healing length £. This suggests that this state can be 
regarded as a topological Wigner crystal. Yet another 
difference is related with the absence of the logarithmic 
divergence of the Coulomb potential ~ \jx leading to 
a very slow decay in the 4fc^ correlation (much slower 
than any power law QI3). Excitations about the qWC 
equilibrium corresponds to the topological equivalent of 
the so-called charge density wave, occurring for K < 1/4 
and characterized by commensurate condition of the fill- 
ing factor 7 = 1/v, where v is an integer. In Fig. ^ we 
illustrate the relevant phases occurring in the system. We 
can observe that it exhibits the same phases as in elec- 
tronic systems, except for the gaseous phase occurring for 
very dilute systems. This stems in the short-range na- 
ture of the topological interactions between the HSs. We 
performed numerical simulations of the Gross-Pitaevskii 
equation to illustrate the existence of both crystal- 
lized and liquid phases. In these simulations, the initial 
configuration was a condensate with a Fermi distribu- 
tion of HSs. Fig. Q depicts the circular polarization 
degree p c = (n+ — n_)/(n+ + which is the most 
suitable quantity to track the dynamics of HSs in BECs. 
In one case (panel a), the regular pattern corresponding 
to Wigner crystal is conserved while HSs are oscillat- 
ing around their equilibrium positions; in the other case 
(panel b), the pattern changes dramatically and the or- 
dering is lost, as it is typical of a liquid-like state. 

Notice that the present analysis is valid in zero tem- 
perature limit with no disorder. Beyond the obvious 
modification of the distribution, it is known that both 
temperature and disorder can melt the crystal state |30j . 
In analogy with the criteria used in one-component plas- 
mas [31], the melting occurs for T = E-^/E^ ~ 170. 
Therefore, we should expect the crystal to be robust for 
HA ^> -Ejnt, -Edis, where E^is is the disorder amplitude. A 
more precise estimation is out of the scope of the present 
work. 

In conclusion, we have investigated the crystalliza- 
tion of a gas of HSs in a spinor Bose-Einstein conden- 
sate. Starting with a kinetic equation to describe the 
coupled dynamics of two-spin species mixture in Bose- 
Einstein condensate, we have shown that the system sus- 
tains unstable oscillations about the Fermi gas equilib- 
rium and spontaneously undergoes crystallization, as a 
consequence of the competition between the interactions 
and the Fermi statistics. For the ordered phase, we com- 
puted the spectrum of collective excitations, obtaining 
two phonon modes (corresponding to the acoustic and 
optical branches). The analysis of the low-energy fluc- 
tuations based on the Luttinger liquid theory revealed 
that the ordered phase indeed corresponds to a Wigner 
crystal, where the Akp correlations are induced by the 
topological effective potential between the solitons. 
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FIG. 3: (color online) Phase diagram of the ordered phase 
described by the Hamiltonian (TTJ. The contours K — 1/2 
(blue line) and K = 1/3 (purple line) delimit the Luttinger 
Liquid (LL), the topological quasi-Wigner crystal (qWC) and 
the mixture (M) phases. Inset: Cut the Luttinger parameter 
K along 7 + = 7_ = 7 for single (dashed) and two-species 
(solid) crystal. The small dashed region represents a zoom- 
out of the gas phase of Fig. |l]). 
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FIG. 4: (color online) Circular polarization degree of the con- 
densate as a function of time showing (a) a stable Wigner 
crystal (7+ = 7_ « 0.25), and (b) "melting" of such crystal 
into a liquid- like state, (7+ = 7_ ~ 0.5). In both situations, 
we took | A| = 0.1. 
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